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ABSTRACT 

We present the BaLROG (Bars in Low Redshift Optical Galaxies) sample of 16 morpholog¬ 
ically distinct barred spirals to characterise observationally the influence of bars on nearby 
galaxies. Each galaxy is a mosaic of several pointings observed with the IFU spectrograph 
SAURON leading to a tenfold sharper spatial resolution 100 pc) compared to ongoing IFU 
surveys. In this paper we focus on the kinematic properties. We calculate the bar strength 
from classical torque analysis using 3.6 /im Spitzer (S^G) images, but also develop a new 
method based solely on the kinematics. A correlation between the two measurements is found 
and backed up by N-body simulations, verifying the measurement of Q^. We And that bar 
strengths from ionised gas kinematics are ~2.5 larger than those measured from stellar kine¬ 
matics and that stronger bars have enhanced influence on inner kinematic features. We detect 
that stellar angular momentum "dips" at 0.2±0.1 bar lengths and half of our sample exhibits an 
anti-correlation of hs - stellar velocity (v/a) in these central parts. An increased flattening of 
the stellar a gradient with increasing bar strength supports the notion of bar-induced orbit mix¬ 
ing. These measurements set important constraints on the spatial scales, namely an increasing 
influence in the central regions (0.1-0.5 bar lengths), revealed by kinematic signatures due to 
bar-driven secular evolution in present day galaxies. 

Key words: galaxies: kinematics and dynamics, galaxies: evolution, galaxies: formation, 
galaxies: structure, galaxies: bulges, techniques: spectroscopic 


1 INTRODUCTION 

The elongated shapes of bars in nearby galaxies were already iden¬ 
tified by Hubble in the mid-30’s ( |Hubble|[l936| l. By the 90s, re¬ 
ported bar fractions reached values of ~30% in local galaxies, each 
in SBs and SABs ( |de Vaucouleurs et 'n^|1991| ). With the arrival 
of extensive imaging surveys, several works reported that at least 
2/3 of the local disc galaxies exhibit bars (e.g., [Knapen, Shlos-| 
[man & Peletier| [20001 |Eskridge et al.|[2000{ |Sheth et al.||2008| l, 
with an increasing fraction found in the infrared. Our own Galaxy 
was also found to host a bar (e.g. [Blitz & Spergel|[T991| t, con¬ 
firmed by large spectroscopic surveys detecting prominent cylindri- 


cal rotation (e.g. BRAVA [Howard et al.|2009[ ARGOS [Ness et al. 

2013 

1 . Bars are ubiquitous in the local Universe (e.g. [Eskridge et al. 

2000 

[Knapen, Shlosman & Peletier[2000[[Whyte et al.[2002[ Mari- 

nova 

& Jogeej2007| jMenendez-Delmestre et al.|2007|[Barazza, Jo- 

gee & Marinova[[2008| [Aguerri, Mendez-Abreu & Corsini[[2009| 

Mendez-Abreu, Sanchez-Jans sen & Aguerri[[2010| [Masters et al.[ 
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[201 1[ [Mendez- Abreu et al.||2014[ [Cisternas et al.[[2014[ l and have 


even been found at higher redshifts (z) (e.g., Abraham et al.|199^ 


Elmegreen, Elmegreen & Hirst 2004| [Jogee et al.|2004|[Simmons[ 

~l.|2014i r but their high-^ fraction seems to be smaller than in 
the local Universe (e.g., [Sheth et al.|2008][Nair & Abraham|2010[ |. 
But how do bars form and how do they infiuence the evolution of a 
galaxy? 


In the early Universe, galaxies were closer together and their 
evolution was dominated by interactions. Nowadays, due to the 
accelerating expansion of the Universe, distances between galax¬ 
ies are increasing. Therefore, the internal evolution experienced by 
the galaxy in isolation - usually termed secular evolution - is in¬ 
creasingly important. Theoretical studies propose that bars have a 
great influence on this internal evolution (see, e.g., [Kormendy &[ 
[Kennicutt|[2004[ for a review). N-body simulations have helped to 
understand barred galaxies, their orbital structures and their influ¬ 
ence on different gala xy properties (e.g. [Combes & Sanders|198l1 
Athanassoulal 1992a|b[ [Debattista & Sellwood|2000[ [Athanassoula 

2003 [ [Bureau & Athanassoul^[2005[ [Martinez- Valpuesta, Shlos- 
man & Heller|2006[ [Minchev & Faniaey|2010[[Wang et al.|2012[ 
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2 Seidel et al. 


[Athanassoula, Machado & Rodionov|2013^ . Because of their sig- 
nificant departure from axisymmetry and the associated torques, 
bars are likely to play a key role in disc galaxy evolution via a 
number of processes: (1) redistribution of the dissipative compo¬ 
nent (as witnessed by e.g., central mass concentration, |Sakamoto| 
|et al.|1999[|Knapen et al.|1995| or flattening of chemical abundance 
gradients, Martin & Roy|1994] l; (2) triggered star formation, associ¬ 
ated with bar-driven crowding and shocks (e.g. spiral arms, rings) or 
central gas accumulation (nuclear starbursts and AGN, e.g. |Coelho| 
& Gadotti|2011|l ; (3) redi stribution of the stellar component (e.g.. 


Gadotti & dos Anjos|2001| l; (4) heating of the stellar component and 


build-up of bulg e structures (e.g., [Bureau & Freeman|199'^|Chung| 
|& Bureau|2004[[Fathi & Peletier|2003|l! 

Many observational studies have investigated barred galaxies 
in profound detail as well as statistically and provided insights 
on their photometry, morphology, kinematics and stellar popula- 


tions (e.g., Sheth et al.||2005| [Gadotti & de Souza||2006| 

|Lau-| 

rikainen et al.[[2007| [Buta et al.[[20101 [Laurikainen et al.f 

2011 

Perez & Sanchez-Blazquez[[2011| [Sanchez-Blazquez et al.[ 

2011 

2014| [Barrera-Ballesteros et al.[[2014| [Kim et al.[[2014[). Because 


of the relative ease with which gas can be observed, the gas dynam¬ 
ics of barred galaxies are fairly well known, strongly constraining 
the models available (hydro, SPH, sticky particles). In contrast, the 
stellar component has not been studied extensively yet, so that the 
dynamical theories (potential and orbital structure) remain poorly 
tested. In particular, our study aims to perform a thorough and direct 
test of the increasing effect on the stellar kinematics with increasing 
bar strength. The consistency of bar-driven bulge formation scenar¬ 
ios (through gaseous inflow, recurring bar formation/destruction, 
recurrent buckling events or other similar mechanisms) also still 
needs to be understood in more detail. 

One of the most important parameters in the characterisation 
of bars (together with properties such as bar length or pattern speed) 
is its "strength". Since our aim is to quantify the influence of the bar 
on the host galaxy, we start our study with a detailed determination 
of this parameter. Analysing our sample in this way, we hope to And 
indications of common features that will be more pronounced in 
more strongly barred systems. Numerous attempts have been made 
to deflne a bar strength parameter in the past, such as the bar axis ra¬ 
tio (e.g., |Martin|1995[[Martinet & Friedli|1997| t. |Aguerri, Beckman| 


|& Prietoj i 1998| l rel^d the amplitude of m = 2 and m = 0 com¬ 
ponents as a measure of bar strength. In 1981, [Combes & Sanders] 
( |1981j ) suggested a measure of bar strength, based on the maximum 
of the bar induced tangential force, normalised to the axisymmetric 
radial force held, denoted as , and widely used in observational 
studies (e.g .,|Buta & Block|2001[|Laurikainen & Saloj200^[Block] 
[et al.|2004}[Buta et al.|2005[[Salo et al.|2010[|. To cal culate in 

this study, we use the polar method {[Salo et al.|1999| [Laurikainen] 

[& Salo[2002] . 


Numerical simulations (e.g.. 

Athanassoula|2003[ [2005||Bour- 

naud. Combes & Semelin[[2005| 

Kim et al.[[2012| [Kim & Stone 

2012|[Lokas et al.|[2014[l often use comparable measurements for 


bar strength and And correlations of increased bar strength with in¬ 
creased angular momentum exchange, higher masses and velocity 
dispersions, but also investigate the influence of magnetic flelds and 


dust lanes (e.g., [Peeples & Martini[2006] l. In this paper we introduce 
a new method of measuring the bar strength using stellar velocity 
maps, based on the recent kinematic decomposition developed by 
[Maciejewski, Emsellem & Krajnovic| ( [2012[ ). This will provide an 
independent view and check of the assumptions made when com¬ 
puting using only the photometry. 

In order to quantify the influence of bars on the properties 


of their host galaxies, we present in this paper integral-fleld spec¬ 
troscopic observations of a sample of 16 nearby barred galaxies. 
The panoramic fleld of view (FoV), which we combined to large 
mosaics for each galaxy, allows us to investigate these galaxies in 
unique detail and probe the influence of bars from the inner regions 
up until co-rotation. This is the first of a series of papers trying to 
unravel the influence of bars in detail combining this high spatial 
resolution dataset with recent simulations and novel analysis tech¬ 
niques. In this paper we will focus on the kinematics and establish 
a yardstick to measure the strength of the bars in order to test their 
influence on other parameters. In a forthcoming paper we will con¬ 
tinue with a stellar population analysis (BaLROG II, Seidel et al. in 
prep.). 

The paper is organised as follows. In section|^we explain the 
motivation of our sample, the observations and data reduction. Sec¬ 
tion summarises the methods employed, particularly explaining 
our calculation of the kinematic torque. We present our results in 
Section and discuss their implications in Section Our main 
conclusions are summarised in Section Throughout the paper 
the cosmological parameters used are Ho = 67kms“^ s“^ Mpc“^, 
Oyv = 0.7 and Orn = 0.3 ( [Planck Collaboration et al.[2014[ l . 


2 DATA 


2.1 Sample selection 

The parent sample from which our target galaxie s are drawn is the 


S^G survey of nearby galaxies ( [Sheth et ah 


20101. We restricted our 


initial choice to barred galaxies with inclinations below 70° and 
brighter than = —18.0 mag to ensure high quality data. As 
we took our sample from the S^G survey, we were naturally con¬ 
strained to those galaxies with cz ^ 3000 km s“^ so that important 
spectral features (H/3, Mgb) remain within the wavelength range 
probed by SAURON. This instrument is mounted on the William 
Herschel Telescope (WHT) in La Palma at the Observatorio del 
Roque de los Muchachos. Therefore, galaxies with sky declinations 
between —2° ^ 5 ^ 60° were chosen to achieve optimal visi¬ 
bility. The entire exercise resulted in a large number of galaxies 
(~ 100), most of which were located in the vicinity of or within 
the Virgo Cluster. Full 2D spectroscopic analysis of a large sample 
was beyond our capabilities in terms of observing time. We thus 
carefully inspected different sets of numerical simulations and im¬ 
ages of the S^G survey and selected those galaxies with prominent 
bars in different apparent stages of their evolution and with dif¬ 
ferent morphologies. In addition, we selected both early-type and 
late-type galaxies and those with inclinations below 60° to reduce 
uncertainties (e.g. in the determination). 

Our target sample consists of 16 galaxies (see Fig.[^, a num¬ 
ber that provides a reasonable representation of different types of 
bars. The number of SAURON pointings greatly exceeds this num¬ 
ber: the dataset for each galaxy is a mosaic of several pointings (up 
to seven individual IFU pointings) allowing us to reach the spatial 
detail that we aimed for while also covering the bars out to the be¬ 
ginning of their surrounding discs. The limitations of the size of 
our sample are obvious. However, while large integral held surveys 
such as CALIFA, SAMI or MaNGA provide large enough samples 
for statistics, they lack the detailed spatial sampling provided by 
this work (e.g. we sample at typically 100 pc, even maintained with 
our Voronoi-binning (within the bar region) which is in most cases 
at least a factor 10 better than the larger surveys). Tablegives the 
entire list of observed targets and basic properties. 
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The BaLROG project -1. Bar influence on kinematics 3 



Figure 1. Spitzer images drawn from the S'^G survey of nearby galaxies jSheth et al.|2010) of our sample observed with SAURON. The number of IFU 
pointings is indicated in the left lower part for each galaxy. The final SAURON mosaic, composed of individual pointings of 30x40 arcsec FoV, is shown on 
top of the images. North is up and East is towards the left in all cases. The extensions of the S^G fields shown here (squares with indicated mosaic on top) are 
160 X 160 arcsec, except for NGC 2893 where it is only lOOx 100 arcsec. 


2.2 Observations 

The observations were carried out in four consecutive runs in March 
2012, January 2013, April 2013 and January 2014 at the WHT 
in La Palma with the SAURON integral field unit { [Bacon et al.| 
|2001j l. In the employed low-resolution mode, this instrument has 
a field of view of 33" x 41", spatial sampling of 0"94 x 0"94 per 
lenslet (1431 in total) and a spectral resolution of full width at half 
maximum (FWHM) of 3.9 A. Its wavelength coverage ranges from 
4760 A to 5300 A, leading to the above mentioned redshift limita¬ 
tions chosen to include important emission and stellar absorption 
line features. 

We observed up to 7 SAURON pointings per galaxy to build 
a large mosaic. The final maps extend along the bars up until the 


start of the disc, allowing us to probe radial dependencies within 
and outside the bar, while also resolving great spatial detail. This 
strategy was quite costly in time: for the small sample of 16 galaxies 
we invested 54 pointings in total, each of 1-2 hours depending on 
the galaxy’s surface brightness. Table [^summarises the number of 
pointings and the total exposure times for each galaxy. Figure 
shows the final extent of the mosaic overlaid on top of the S^G 
images. 

Apart from the large pointing offsets, we introduced small 
dithers within each pointing of typically 1 to 2". This helps us to 
account for a couple of bad columns in the CCD and to improve our 
sampling. The orientation of the field of view (FoV) of SAURON 
was such that the 146 sky lenslets always pointed away from the 
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Table 1. Galaxy Sample. - (1) Galaxy name, (2) Hubble type, (3, 4) J2000 coordinates (right ascension, declination), (5) systemic velocity, (6) inclination, (7) 
bar length, (8) position angle, (9) bar position angle, (10) effective radius (rg). Objects forming part of the Virgo cluster are marked with a small v next to their 
name. Notes. - All morphological classifications, coordinates and Vgys are from the NASA Extragalactic Database (NED), http://ned.ipac.caltech.edu/ All 
inclinations, PAs and effective radii of the galaxy (rg) are from the S'^G P4 (Salo et al. 2015, accepted to ApJS) while bar lengths and bar PAs are determined 
by Herrera-Endoqui et al. (2015, submitted) (here the ones by visual inspection). 


Galaxy 

(1) 

Hubble 

Type 

(2) 

RA 

(hh mm ss.ss) 
(3) 

Dec 

(dd mm ss.ss) 
(4) 

Vsys 

(kms“^) 

(5) 

Inclination 

(deg) 

(6) 

Bar Length 
(arcsec) 

(7) 

PA 

(degrees) 

(8) 

FA-bar 

(degrees) 

(9) 

re 

(arcsec) 

(10) 

NGC 1015 

SBa 

02 38 11.56 

-01 19 07.3 

2628 

30.5 

21.5 

-8.3 

101.4 

19.73 

NGC 2543 

SB(s)b 

08 12 57.92 

+36 15 16.7 

2471 

59.9 

14.9 

37.0 

105.5 

26.56 

NGC 2712 

SBb 

08 59 30.47 

+44 54 50.0 

1815 

60.5 

20.5 

3.6 

22.6 

25.09 

NGC 2859 

(R)SB0+(r) 

09 24 18.53 

+34 30 48.6 

1687 

37.2 

34.4 

1.8 

169.6 

22.41 

NGC 2893 

SBO-a 

09 30 16.96 

+29 32 23.9 

1703 

17.4 

12.0 

88.1 

164.1 

4.58 

NGC 2962 

SBO-a 

09 40 53.93 

+05 09 56.9 

1960 

49.0 

30.7 

8.3 

172.7 

20.50 

NGC 3485 

SBb 

11 00 02.38 

+14 50 29.7 

1436 

20.4 

21.0 

-64.6 

40.6 

26.38 

NGC 3504 

SBab 

11 03 11.21 

+27 58 21.0 

1539 

12.8 

37.1 

-41.7 

148.3 

11.13 

NGC 4245 

SB0/a(r) 

12 17 36.77 

+29 36 28.8 

886 

33.3 

36.3 

0.5 

131.0 

23.52 

NGC 4262^ 

SBO-(s) 

12 19 30.57 

+14 52 39.6 

1359 

24.5 

13.4 

-6.0 

26.5 

5.99 

NGC 4267^ 

SBO- 

12 19 45.24 

+12 47 53.8 

983 

11.9 

16.9 

-27.5 

34.0 

21.07 

NGC 4394^ 

(R)SB(r)b 

12 25 55.53 

+18 12 50.6 

922 

30.4 

41.4 

-57.6 

143.4 

36.79 

NGC 4643 

SB0/a(rs) 

12 43 20.14 

+01 58 41.8 

1330 

36.8 

49.9 

56.0 

133.3 

24.22 

NGC 5350 

SBbc 

13 53 21.63 

+40 21 50.2 

2321 

50.3 

15.2 

7.9 

120.8 

28.06 

NGC 5375 

SBab 

13 56 56.00 

+29 09 51.7 

2386 

29.8 

27.2 

-9.4 

171.1 

24.35 

NGC 5701 

(R)SB0/a(rs) 

14 39 11.08 

+05 21 48.5 

1505 

15.2 

39.0 

52.0 

174.9 

25.97 


Table 2. Summary of the observations: (1) NGC number, (2) Run number, 
(3) Pointing number, (4) Total exposure time, in seconds. 


Galaxy Run # P Tgxp 

(1) (2) (3) (4) 


NGC 1015 2 

NGC 2543 1,2 

NGC 2712 4 

NGC 2859 1 

NGC 2893 3 

NGC 2962 4 

NGC 3485 4 

NGC 3504 3 

NGC 4245 2 

NGC 4262 1 

NGC 4267 1 

NGC 4394 2 

NGC 4643 3 

NGC 5350 2 

NGC 5375 4 

NGC 5701 1 


2 12 X 1800 

3 10x1800 

4 8 X 1800 

4 16 X 1800 

1 8 X 1800 

3 9 X 1800 

2 4 X 1800 

4 12 X 1800 

4 16 X 1800 

2 8 X 1800 

3 12 X 1800 

6 24 X 1800 

7 23 X 1800 

2 8 X 1800 

4 16 X 1800 

3 12 X 1800 


galaxy’s centre. They are 1.9 ' from the main FoV and thus ensure 
a simultaneous sky exposure during the object exposure. 

We took a calibration frame using a Neon lamp before and 
after each science frame. Skyflats were taken at dusk and dawn, 
as well as continuum lamp exposures with a Tungsten lamp. For 
the flux calibration we observed several spectrophotometric stan¬ 
dard stars. For further spectral calibration, we also observed a 
broad rang e of stars with different spectr al types from the MILES 
databas^ ( [Sanchez-Blazquez et al.|200^ . 


^ http://miles.iac.es 


2.3 Data reduction 

The reduction was performed with the available SAURON pipeline 
XSauron described in detail in [Bacon et al.| ( |2001| ). The prepro¬ 
cessing of raw frames includes overscan and bias subtraction. The 
evaluation of dark frames showed that the dark current is negligi¬ 
ble: less than 1 e“pixel“^h“^. All frames were preprocessed in this 
same way. After that, a model mask was created to extract the spec¬ 
tra. This mask builds a table indicating corresponding positions by 
relating the pixels on the CCD to their associated wavelengths and 
lenslets. The outcome is a set of three-dimensional data cubes (a, S, 
A). Wavelength calibration was achieved with the arc (neon) lamp 
exposures. A cross-correlation function between the neon frames 
taken before and after the science exposure and the one of the ex¬ 
traction mask deflnes potential slight offsets between the science 
frame and the mask. This analysis is based on 11 emission lines 
which can be seen in the wavelength range of SAURON. 

The flat-flelds were created with a combination of twilight and 
continuum lamp (tungsten) exposures. The former calibrates the 
spatial component, while the latter is responsible for the spectral 
coordinates. We used for each run a representative twilight and con¬ 
tinuum flat exposure investigating counts and distributions of all flat 
exposures. On a case-by-case basis, we also chose night-dependent 
flats, but for the vast majority and thanks to our bright objects, the 
former method proved to work well. Cosmic rays were removed 
before the sky subtraction, where the median of the 146 dedicate 
sky lenslet values was computed and subtracted from the science 
frame spectra. Flux calibration was done using the spectrophoto¬ 
metric standards. Their flux correction curve was extracted com¬ 
paring the observed curve with a reference spectrum. The resulting 
correction curve was used to calibrate all science frames. The merg¬ 
ing and mosaicking of the individual data cubes was achieved with 
the XSAURON software using the integrated intensity contours in 
comparison with those of a ^-band SDSS image. The entire mo¬ 
saic was constructed with the obtained offsets and scalings between 
each image. 
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2.4 S^G data 


We complement ou r SAURON mosai cs with Spitzer 3.6/im im¬ 
ages from the S^G ( |Sheth et al.|2010 i. As the Spitzer images are 
very deep, the outer isophotes are typically close to or beyond the 
3.6 = 25.5 mag/arcsec^ in the AB magnitude system, the position 
angles (PA) and ellipticities {i) are taken from the S4G pipeline 4 
(Salo et al. 2015, accepted to ApJ). Global galaxy parameters such 
as the effective radii of the galaxies (re) are from Munoz-Mateos 
et al. (2015, submitted) and the barlength measurements are from 
Herrera-Endoqui et al. (2015, submitted). In addition, we used the 
Spitzer images to compute bar strengths for our sample as de¬ 
scribed in Section |3.3.1| The measurements for the complete 
S^G are given in Diaz-Garcia et al. (2015, submitted). 


3 METHODS 


A detailed analysis of the stellar and gas kinematics requires a mini¬ 
mum signal-to-noise ratio (S/N) (e.g., |van der Marel & Franx| 1993|l. 
We adopted the Voronoi binning scheme of Cappellari & Copin| 
( |2003| ) and applied it to our data to reach a minimum S/N?^40 
per pixel for all galaxies. The central spectra remained unbinned 
in all cases and exceeded this S/N level (e.g. S/N>100) while in 
the barred regions, we reach a typical S/N of . While a S/N of 40 
ensures high-quality spectra for the extraction of the mean stellar 
velocity, velocity dispersion as well as Gauss-Hermite moments ha 
and h 4 , it is also low enough to preserve the spatial substructures 
in the galaxies, as seen in the resulting maps (see Appendix [^. 
Before we binned, we also ensured that we would not contaminate 
our measurements by poor quality spaxels. Therefore we excluded 
those spaxels with a S/N below 3 and then limited the data to an 
isophote with at least this average S/N level. The resulting exten¬ 
sions of the maps are hence due to the combined mosaic and this 
additional S/N minimum threshold. 


3.1 Stellar kinematics 


We extracted the stellar kinematics using the pPXF - penalized 
pixel fitting - code developed by |Cappellari & Fmsellem| ( |2004| ). 
The routine fits each galaxy spectrum with a combination of tem¬ 
plate spectra from a given library. Here we use a subset of Medium- 
resolution Isaac Newton Telescope library of empirical spectra 
(MIFFS; [Sanchez-Blazquez et al.||20 06|l s ingle stellar population 
(SSP) model spectra ( Vazdekis et al.|20]^ with a a range of ages 
and metallicities of 0.1 Gyr to 17.8 Gyr and -0.40<[Z/HI <-h0.22, 
respectively. Their mean resolution is of FWHM = 2.51 A |Falc6n 
[Barroso et akpOllT l and before the fitting process, we matched the 
spectral resolution of the models to that of our data. Throughout 
this work we assume a Kroupa initial mass function (IMF, |Kroupa| 
|2001| l. The result of pPXF is a line-of-sight velocity distribution 
(LOSVD) described by a Gauss-Hermite parametrisation { |Gerhard| 
|1993[ |van der Marel & Franx|199^ allowing the measurement of 
the velocity (I/), velocity dispersion (cr) and higher order Gauss- 
Hermite moments (ha and h 4 ). From the stellar velocity and veloc¬ 
ity dispersion, we also calculated the value of the specific stellar 
angular momentum Ar ( [Fmsellem et al.|[2007| ), radially but also 
integrated within Ire. Using the code developed by |Maciejewski,| 
|Fmsellem & Krajnovic| ( |2012] t, we furthermore obtained the radial 
and tangential velocities Vr and Vt for a subset of our sample (see 
Section |J.3.2| for more details). 



Figure 2. Left panel: Gravitational torque map of NGC 2543 derived from 
Spitzer image using only even Fourier components. In the top of the panel 
a colour bar shows the maximum and minimum Qq; Bar length and ellip- 
ticity are traced with a black solid line; The inner dotted circle corresponds 
to the S'^G Pipeline 4 re of the bulge. Right panel: Fourier smoothed density 
with the axisymmetric component (m = 0) subtracted; Contours of equal 
Qg are overplotted in white; The dotted lines indicate the regions where the 
tangential forces change sign. 


3.2 Emission lines 

We used the Gas AND Absorption Line Fitting (GANDALF) pack¬ 
age by |Sarzi et al.| ( |2006T ) and |Falc6n-Barroso et al.| ( [2006] ) to deter¬ 
mine the ionised-gas distribution and kinematics. The code treats 
the emission lines as additional Gaussian templates on top of the 
stellar continuum and iteratively searches for the best match of 
their velocities and velocity dispersions. The SAURON wavelength 
range allows us to measure the emission line of H/3AA4861 and 
the doublets [Oiii]AA4959, 5007 and [Ni]AA5200, 5202 A. We tied 
spectral lines kinematically to the [Olll] doublet to lower the num¬ 
ber of free parameters given to GANDALF. We checked that leaving 
them free resulted in a consistent outcome. Furthermore, we im¬ 
posed known relative fiux relations to constrain the freedom of the 
doublet lines during the fitting process, namely F([Olll] 4959 ) = 
0.350 • F([Oiii]5oo7). 

We will use the results from GANDALF to clean the spectra 
and produce emission-line-free datacubes for our stellar population 
analysis, to be presented in our Paper II. 


3.3 Bar strength measurements 

Bar strengths have been measured in many different ways (see in¬ 
troduction for details). For our analysis we will concentrate on the 
following two methods: (1) the photometric torque taking ad¬ 
vantage of the S^G data and (2) a new measurement based on the 
stellar velocity maps which does not include strong model assump¬ 
tions (Qkin). 


3.3.1 Photometric torque using 3.6 pm Spitzer imaging (Q^) 


We calculated the gravitational potential of our galaxies from the 
3.6 micron images, using the NIRQB code ( [Laurikainen & SaIo| 
|2002| t based on the polar method developed in Salo et al. \1999\ . 
Before applying the Fourier transformation, the Spitzer images are 
rectified to face-on. Then, the even Fourier components (up to 20) 
of the surface density / (r, f) are calculated within a polar grid. The 
gravitational potential 0) is then inferred from the smoothed 

surface densities by applying a fast Fourier transformation in the 
azimuthal direction in combination with a direct summation over 
radial and vertical directions. We use a polar grid with 128 bins in 
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the azimuthal direction, which determines an angle step-size for the 
azimuthal Fourier transform of — 2.8°. 

The calculation of the potential is based on the following as¬ 
sumptions: 

(i) The mass-to-light ratio is constant. 

(ii) The disc vertical scale height hz is constant. 

(iii) The disc has an exponential vertical density distribution: 

Pz{z) = ^^exp{-\z/h:,\). (1) 

(iv) The vertical scale height of the disc scales with the disc size 
as hz =0.1 rk 2 o» where rk 2 o the 2MASS { [Skrutskie et al.|20()6| 
iT—band surface brightness isophote of 20 magarcsec”^. 

Tangential forces {FT{r^ and radial forces 

(FR(r, 0)/^r) are obtained via integration. Non- 

axisymmetric forces in the galaxy are characterised by the ratio of 
the tangential force to the mean axisymmetric radial force field: 

QG(r,0) = FT(r,0)/(FR(r)), (2) 

where (FR(r)) is the azimuthally averaged radial force at a radial 
distance r. Qg(^, 0) values are used to construct the gravitational 
torque maps of our galaxies (see example in Fig. [^. Typically, 
barred galaxies show a well-defined four-quadrant Qg map, resem¬ 
bling a butterfly pattern, which is roughly symmetric with respect to 
the bar major axis. We take the even Fourier components uniquely 
(focusing on bi-symmetric structures) and we symmetrise our maps, 
reducing in this way the impact of sharp density clumps. 

Based on the torque maps, and given a certain radial dis¬ 
tance r and quadrant q, one can identify a maximum QT(r)^ = 
max(QG(r, 0)^). We calculate the radial profile of the relative 
strength of the non-axisymmetric perturbations throughout the 
galaxy, QT(r), taking the mean of these four maxima. For addi¬ 
tional information about the method, see |Salo et al.| ( |20l6] ). Finally, 
the gravitational torque parameter (Q^) corresponds to the maxi¬ 
mum value of Qt at the bar region. 

The main source of uncertainty (~15%) is the poorly known 
vertical thickness: to account for this we have used different disc 
thicknesses in the calculation of the gravitational field. A small sys¬ 
tematic error is produced by the omission of the dark halo contri¬ 
bution on the radial forces, but this is likely to be smaller than that 
associated with the vertical thickness (Diaz-Garcia et al., submit¬ 
ted.) 


3.3.2 Kinematic torque (Qi^in) 

In order to perform a model-independent measurement and to test 
the torque measure of we developed a new method solely using 
the kinematics, resulting in a new parameter which we call the kine¬ 
matic torque Qkin- The basis of this analysis is the stellar velocity 
field. Using this map, we extracted the radial and tangential veloc¬ 
ities following [Maciejewski, Emsellem & Krajno^^ ( |2012^ , using 
their equations 9 and 10. This method is based on assuming a thin 
disc geometry to obtain the two velocity components in the equa¬ 
torial plane. Further assumptions in deriving these two quantities 
are: 

(i) A steady state bar, hence not in buckling phases or alike. 

(ii) A symmetric bar with respect to its major axis. 

(iii) A thin galaxy disc resulting in only two velocity compo¬ 
nents. 



Vradial 


NGC 



Figure 3. Radial and tangential velocities for NGC 2543. Upper left: radial 
velocity field, upper right: tangential velocity field; in both of them, the 
continuous circle indicates the bar radius and the dotted circle the radius 
where we measured the kinematic torque (Qkin)- Lower left: radial velocity 
along the radius; lower right: radial distribution of the tangential velocity; in 
both: grey crosses represent individual measurements, red points show the 
obtained modulation (not a fit to the gray points - see text for details) and 
the blue a smoothed fit to it. The dashed lines indicate the maximum found 
in the radial velocity and corresponding to the same radius shown for the 
tangential velocity (see text for details). 


As the buckling phase is only a brief evolutionary stage of the 
bar (e.g.|Shen & Sellwood|2004||Martinez-Valpuesta, Shlosman & | 

|Heller|2006> nt is much more likely to observe bars in their steady 

state. Under the presence of spiral structure, the bar will lose its 
symmetry with respect to the major axis, but only at the very edges 
of the bar. An aspect to consider is that prominent bulges will break 
the third assumption in the inner regions of the bar. Following these 
limitations, one can conclude that the most reliable region is within 
the middle of the radial extension of the bar, where we expect to 
measure the strongest radial velocities and corresponding torques, 
as outlined below. In addition to these assumptions, the technique 
developed by |Maciejewski, Emsellem & Krajnovic| ( |2012| ) requires 
the knowledge of the systemic velocity, galaxy inclination, position 
angle of the line-of-nodes, the bar position angle and bar length. We 
estimated those from the literature and close inspection of our own 
datasets. Furthermore, we rectified our maps to face-on in order to 
apply the technique. 

Using the extracted radial and tangential velocities, we defined 
a new parameter (Qkin) that measures the torque directly from the 
observed kinematics: 


_ max(urad(R)) 
(abs(utan,R)) ’ 


(3) 


where we first find the radial position of the maximum value of the 
radial velocity (Urad), and then determine the corresponding tan¬ 
gential velocity (utan) as the mean value in a ring around this ra¬ 
dius. This relation is constructed analogous to the calculation of 
based on the fact that Urad is proportional to Ft/Fr x Urot and 
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t'tan roughly equal to v^ot- Therefore Urad/^tan is expected to be 
proportional to Ft /Fr (note that the ratios do not have to be equal, 
but only proportional). Figureshows the radial (left) and tangen¬ 
tial (right) velocity maps (top) and radial (bottom) distribution for 
NGC 2543 as an example. The position of the maximal radial ve¬ 
locity is found by evaluating the radial velocity field in rings. We 
expect a certain velocity modulation when tracing a circle through 
the four quadrants, i.e., combination of sine and cosine curve when 
tracing the radial velocity in a ring. This additional aspect helps 
us to detect and correct for outliers, i.e., unreal peaks or drops of 
extremely high or low values, which appear more often in the kine¬ 
matic data due to higher noise levels. Hence, we avoid to simply 
measure the maximum which would lead to an incorrect result. A 
smooth version is then obtained by fitting a polynomial. In the bot¬ 
tom panels, we show the curves obtained when measuring the am¬ 
plitudes of the modulation (red and smoothed fit in blue) compared 
to the individual data points (grey). It is obvious that there is a sig¬ 
nificant scatter among the individual points, but nevertheless a clear 
maximum can be distinguished in the radial velocity profile which 
is well captured with the modulation. 

To further constrain the measurement, we evaluated Qkin 
within the bar region as determined from the S^G images. This is to 
avoid choosing areas where high values appear, either due to higher 
noise levels towards lower surface brightness areas, or due to spiral 
arms. In the example shown in the figure, the strength of the spiral 
arms can clearly be seen: in a central arcsec radius, we detect 
the signature of the bar, but further out, the field does not become 
flat but shows other maxima and minima due to the torques exerted 
by the spiral arms. Those strong, outer values detected in the radial 
velocity - and thus the measured torque - is not due to the bar but 
to the spiral arms in this galaxy. Similar enhancements can also be 
seen in Fig.[^for the computation of . 

As the value of Qkin depends on the input parameters to de¬ 
termine radial and tangential velocities, we chose to determine its 
uncertainty via a set of Monte-Carlo simulations. For each reali¬ 
sation, we chose a random combination of initial values of the in¬ 
clination, line of nodes position angle, bar position angle and bar 
length, all within their uncertainties. As inclination is the most dif¬ 
ficult to determine, we allowed an uncertainty of ±10°, whereas 
we chose ±5° for the other parameters, leading to an overall uncer¬ 
tainty found in Tab. Higher values would result in simply higher 
uncertainties in the measured torques. 

Unfortunately, the determination of Qkin is only possible when 
the kinematic major axis and the bar position angle are neither per¬ 
pendicular nor parallel (at least 5° off, while an angle of 45° would 
be ideal). It is only under those circumstances that the method of 
|Maciejewski, Emsellem & Krajnovic| ( |2012| ) can be applied to com¬ 
pute the required Urad and Vtan- From the 16 galaxies in our sam¬ 
ple, we could only measure the kinematic torque (Qkin) on the fol¬ 
lowing 10 systems: NGC 2543, NGC 2712, NGC 2859, NGC 2962, 
NGC 3504, NGC 4245, NGC 4262, NGC 4394, NGC 5350 and 
NGC 5701. Results are summarised in Tab.|3] 


4 OBSERVED KINEMATIC PROPERTIES 

This section summarises the different parameters extracted from 
the kinematic maps of stellar and ionised gas component. Fig- 
urej^presents two examples of absorption-line stellar velocity maps 
and associated radial profiles along the major and minor axis for 
two galaxies in our sample, NGC 4643 (early-type) and NGC 4394 
(late-type). The complete set of kinematic maps, including ionised 




Stellar velocity, mean PA: 57.77 
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Figure 4. Stellar velocity maps and their associated major axis profile for 
two galaxies. The colour bars on the side each indicate the range of the 
parameter measured. The isophotes shown are derived from the SAURON 
cube reconstructed intensities and are equally spaced in steps of 0.5 mag¬ 
nitudes. The kinematic position angle, based on |Barrera -Ballesteros et al.| 
( pOU) , is given on the top. The dots indicate the bins used for determin¬ 
ing that angle. The black lines indicate the photometric position angle, even 
better seen in the following figures, not demonstrating the method. 


gas kinematics and Gauss-Hermite moments ha and h 4 are collected 
in Appendix [C| Overlaid in all maps, we show the isophotes of the 
surface brightness (in mag/arcsec^ with an arbitrary zero point) re¬ 
constructed from the SAURON datacubes and equally spaced in 
intervals of 0.5 magnitudes. In this section, we concentrate on an 
overview of the general kinematic trends observed in our sample. 
We also present the bar strength measurements from these kine¬ 
matics (Qkin), in comparison with the ones derived from the S^G 
imaging (Q^). 


4.1 Stellar and gas kinematics 

We investigate the orientation of the stellar and gas kinematics, 
comparing them to the bar axis, as well as to features that can be 
linked to bar-driven secular evolution. We use the entire maps as 
well as cuts along different axes to better unravel certain features. 

A first glance at the maps shown in Appendix [C] reveals that 
the overall rotation is not strongly affected by the bar (i.e., the kine¬ 
matic major axis remains almost constant as a function of radius as 
determined using the method by [Barrera-Ballesteros et al.||2014| l, 
implying that the bar has not changed the global rotation pattern 
of the galaxies. We do not detect either large velocity twists in the 
line-of-nodes (a kinematic feature observed in simulations). Only 
NGC 2712 and NGC 4394 show small deviations. The absence of 
this feature in our maps may be due to projection effects or simply 
to the limited FoV, because the twist is often visible further out, such 
as in NGC 936 (e.g., [Maciejewski, Emsellem & Krajnovi(5|2012] l. 

Along the kinematic major axis, we do detect in all cases the 
so-called double-hump rotation curve (local inner maximum fol¬ 
lowed by a slight drop and further rise) predicted by simulations 
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Figure 5. Stellar velocity dispersion maps and their associated profile along 
the bar major axis for the same two galaxies as in Fig|^ The colour bars 
on the side each indicate the range of the parameter measured. The contin¬ 
uous black line shows the bar major axis, determined from the photometry 
of Spitzer images while the dashed line indicates the global photometric 
position angle found previously. 


( [Bureau & Athanassoula|2005|l, so far mainly confirmed in edge-on 
systems (e.g., Chung & Bureau|2004| l. This feature can be seen in 
the stellar velocity maps as an enhanced area of high (low) velocity 
values on both sides of the nucleus, but is obviously more apparent 
in the radial profiles. This double-hump is clearly visible in more 
than 60% of the galaxies in our sample and hints towards the ex¬ 
istence of inner discs or rings. Along the minor kinematic axis, we 
also find a similar distortion in the very central parts, visible as a 
small-scale wiggle in the minor axis rotation profile. This profile is 
normally expected to be flat with a value around zero, but it appears 
to be present for all cases. This small feature might indicate a non¬ 
perfect estimation of the global photometric position angle of the 
galaxies. 

The stellar velocity dispersion maps, and radial profiles, of the 
two example galaxies NGC 4643 and NGC 4394 are presented in 
Fig.|^ The maps show the presence of cr—drops, a kinematic fea¬ 
ture also predicted by simulations of barred galaxies (e.g., |Wozniak| 
jet aLpOO^ . The fraction of galaxies in our sample that show this 
behaviour is 62.5%. The early types in our sample show a larger 
region of overall higher a than the late types, as expected from the 
presence of larger central bulge structures. Peak velocity dispersion 
values range between 100 km s “ ^ and 220 km s “ ^. 

The corresponding kinematic maps (velocity and velocity dis¬ 
persion) for the ionised gas are shown in Fig. for the two exam¬ 
ple galaxies, illustrating their difference to the stellar kinematics. 
The ionised gas maps of the entire sample appear in general less 
regular with a more patchy distribution, reflecting the gas prop¬ 
erties. We do not find gas equally distributed in all the galaxies. 
The absolute values of minimal and maximal rotation are gener¬ 
ally slightly higher for the gas than for the stars in all galaxies. The 
kinematic major axis of the gas velocity field has the same orienta¬ 
tion as the stellar velocity field. Only one galaxy shows a significant 
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Figure 6. Velocity and velocity dispersion maps of the ionised gas compo¬ 
nent for the two examples in the BaLROG sample of galaxies. 


change in orientation, NGC 4262, whereas two others (NGC 2962 
and NGC 5701) show very mild differences only. Previous studies 
suggest that NGC 4262 might have been involved in an interaction 
( jVollmer, Huchtmeier & van Driel||2005j l. The gas velocity fields 
also present the double-hump feature, in some cases significantly 
more pronounced than in the stellar maps (e.g. NGC 3504). This 
confirms that the gas is more susceptible to bar-driven processes 


(e.g., Schwarz||1981[ [Ellison et al.||201l( [Athanassoula, Machado [ 

[& Rodionov[2013[ l~The ionised gas velocity dispersion differs sig¬ 

nificantly from the stellar a. It does not show a central elevation 
(tracing the bulge component) but exhibits slightly higher values 
throughout the area inside the bar isophotes. However, the pattern is 
extremely patchy and we lack enough coverage of the disc for a fair 
comparison between the morphological components. Especially in 
the late-type systems, some regions display higher velocity disper¬ 
sion in the gas than in the stars. These are typically associated with 
spiral arms (e.g. NGC 3504, NGC 4394). 

We will link the investigated features to the strength of the bars 
in Sec. 15.21 


4.2 Gauss-Hermite moments: ha and \va 

In addition to the first and second moment of the LOSVD distribu¬ 
tion (y, cr), we measured the ha and h 4 Gauss-Hermite moments. 
They help to understand the distribution of orbits along the line-of- 
sight and can be used to distinguish dynamically distinct regions 
and thus indicate whether bars influence their formation. 

Mathematically, ha measures the skewness of the LOSVD, i.e., 
wings on either side of the peak deviating from the otherwise Gaus¬ 
sian profile, while h 4 is a measurement of the kurtosis (e.g. [van d^ 
[Marel & Franx[I99'3] l. Within our sample, we find a large variety in 
those maps and just by visual inspection, we cannot identify a sys¬ 
tematic pattern which could be attributed to the bars of their host 
galaxies. 

Figure |7] shows the two Gauss-Hermite moments ha and h 4 
for NGC 4394 and NGC 4643. In the first galaxy, the maps ap¬ 
pear more uniform than in the second one, which does exhibit low- 
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Figure 8. Individual Voronoi-binned values of ha versus the stellar velocity over the stellar velocity dispersion, within IRf^ar’ 0.5 Rbar and 0.1 
Representative error bars for the red and green regions are indicated in the lower left corner of each panel. 


level structure. However, the majority of the maps of ha and h 4 
(see more in the appendix) show very low values (below 0.1) and 
are rather flat throughout the FoV. In several cases a slight anti¬ 
correlation between the ha moment and the stellar velocity can be 
seen (NGC 1015, NGC2959, NGC2962, NGC4245, NGC4262, 
NGC4643 and NGC 5350), while h 4 moments correlate in most 
cases with the velocity dispersion. We investigated the h 4 profiles 
along the bar major axis for kinematic signatures of peanut-shaped 
bulges (e.g., |Debattista et al.|2005| l, but could not detect any clear 
evidence. This is a property of mainly very low inclined galaxies 
with strong peanut shapes seen in the photometry, so our sample 
is not an ideal selection for the detection of this characteristic. So 
far, only a few studies (e.g., [Mendez-Abreu et al.|2()08]|2014j ) have 
confirmed this prediction. 

In the literature, the ha moment is often related to the stellar 


velocity and both correlations and anti-correlations are found (e.g., 
[Bureau et al.|2004| l. For a more robust measure we chose to corre¬ 
late h3 with V/cr, shown in Fig.|^ in order to compensate for differ¬ 
ent masses. However, we checked the relation correlating with the 
stellar velocity alone and did find very similar results. In the figure, 
Voronoi-binned values within the bar length are shown in black, 
within 0.5 times the bar length in green and within 0.1 times the 
bar length in red. There is no clear (anti-)correlation for the full ex¬ 
tent of the bar length (nor for the effective radius, which we tested 
for comparison, but which is not shown here). However, the more 
central the aperture, the stronger the anti-correlation for about 50% 
of our sample, clearly depicted by the red points. The other ~50% 
of our sample do not show any (anti-)correlation at all, but a simple 
spread of hs values around zero. This behaviour is consistent across 
both measures - bar length or the effective radius - with the excep- 
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Figure 9. Individual Voronoi-binned values of 114 versus the stellar velocity 
the red and green regions are indicated in the lower left corner of each panel. 

tion of NGC 4262. The effective radius for this galaxy is signifi¬ 
cantly smaller and captures the inner part only. Hence, as smaller 
apertures decreased the scatter if an anti-correlation was present, 
it stands out better for the effective radius measurement. In many 
of the other galaxies the effective radius is comparable to the bar 
length or at least not less than half its size. Overall, we cannot de¬ 
tect a tendency between late or early types, because both late types 
(e.g. NGC 4394) and early types (e.g. NGC 2859) show the above 
described behaviour, although with a mild bias towards earlier types 
showing stronger correlations. In the figure, we ordered the galaxies 
according to their Hubble type (SBO top left to SBbc bottom right) 
and one can appreciate the larger scatter also amongst the innermost 
(red) points in the bottom row (latest types of our sample). The fact 
that there is no striking difference though might mean that the Hub¬ 
ble type is not the crucial factor, neither the bar, but the presence of 


dispersion, within 1 -Rbar^ 0.5 Rbar and 0.1 Rbar- Representative error bars for 


significant substructures. More than 50% of those with strong cen¬ 
tral anti-correlations have confirmed substructures such as nuclear 
rings, nuclear lenses or nuclear ring or bar lenses (see classifications 
by |Laurikainen et al.|2011| and |Buta et al.|2015| ). 

The influence of bars on building up a central component is 
supported by this hs-V anti-correlation in the centres of about 50% 
of our sample of galaxies. Earlier studies have found a correlation 
as well as anti-correlation between the stellar velocity and ha mo¬ 
ment, depending on the area and type of galaxy sampled. In edge-on 
barred galaxies, [Bureau et al.| ( |2004| ) and [Bureau & Athanassoula| 
( [2005| ) detected an hs-V correlation over the projected bar length, 
expected for a thick bar. In the centres, however, they also found 
an anti-correlation in more than 60% of the galaxies. This can indi¬ 
cate the presence of multiple components with different kinematics. 
Hence a significant number of barred galaxies, not only in edge- 
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Figure 7. Gauss-Hermite maps of the stellar component for the two exam¬ 
ples in the BaLROG sample of galaxies. 


on systems but also in our sample of different inclinations, show 
the presence of cold and dense (quasi-)axisymmetric central stel¬ 
lar discs. This supports the scenario of the bar driving gas towards 
the centre and nourishing star formation, resulting in this additional 
central component. The coincidence of a steep central light profile 


and star-forming ionised gas discs in these same regions (e.g. Bu- 


reau & Freeman|19^ supports this theory further. [Falcon-Barroso 
et al.| ( 2006| l also found a link between more intense star forming 


regions and lower gas dispersion values. However, due to our lim¬ 
ited wavelength range, we could not account consistently for dust 
correction which could infiuence the values of ha, being sensitive 
to dust. Dust is particularly important in late-type galaxies, but the 
majority of our sample are early-types and no obvious dust lanes 
are detected in the later types. 

Figurej^ shows similar measurements for The galaxies are 
again ordered by Hubble type. This parameter measures the sym¬ 
metric deviation from a Gaussian profile, indicating a velocity dis¬ 
tribution which is less (or more) peaked (negative or positive val¬ 
ues). It is expected to correlate with the velocity dispersion ( |van| 
|der Marel & Franx|1993| l, so we chose to plot it against a. Over¬ 
all, but especially for the early types, we observe elevated values in 
the central regions, most probably associated with the bulge. The 
presence of higher values together with an occasional cr—drop 
hints at the presence of components with more recent star forma¬ 
tion within a classical elliptical-like bulge component { |Wozniak| 
|et aL|2003{ [Bureau & Athanassoula|20()5] l. In Fig. [9] this behaviour 
stands out even more clearly: a higher velocity dispersion is found 
in the centres where the red points (corresponding to most central 
values, i.e., within 0.1 bar lengths) also show elevated values. 
Later-type galaxies are clearly different. Not only does the velocity 
dispersion cover a larger range of values, also the h 4 values are not 
particularly high in the centre or even go significantly below zero 
(e.g. NGC2712, NGC3504, NGC3485, NGC5350, NGC5875). 
The bulge in the later-type galaxies is significantly less pronounced, 
therefore the mixture of more components could cause this spread 
of h 4 values. The strongest central concentration of elevated h 4 
values, at an almost constant a and h 4 , are found in NGC 1015, 


NGC2859, NGC 2962, NGC 4262 and NGC 4267. These are not 
the galaxies showing the ha-velocity central anti-correlation, but 
those to have a prominent bulge component. NGC 5701 also has a 
large bulge, but contains confirmed nuclear spiral structure ( |Erwin| 
|2004| ), which could contribute to the higher spread in h 4 . 

In conclusion, the analysis of Gauss-Hermite moments sug¬ 
gests that the centers of barred galaxies (within at least half a bar 
length, and even more obvious within 0.1 bar lengths) host dynam¬ 
ically distinct components. These could have been altered by bar- 
driven evolution. 


4.3 Angular momentum: Ar 


We calculate Ar as a measure of rotational versus pressure sup¬ 
port following the prescription given in |Emsellem et aL|p007| ), as 
bars are meant to work as engines redistributing angular momen¬ 
tum amongst the different components of a galaxy. This parameter 
is based on the first two stellar velocity moments and the corre¬ 
sponding fiux and is defined as 




(4) 


for two-dimensional spectroscopy, where Fi denotes the fiux, Ri 
the circular radius, Vi the velocity and ai the velocity dispersion of 
the zth spatial bin (going to Np bins). As outlined in |Emsellem et al.| 
( |2011| ), it improves the characterisation of the dynamical state of a 
galaxy compared to the simple measure ofV/a. It shows a clear 
difference, especially for non-regular rotators with irregularities in 
their velocity fields, whilst at the same time being correlated to the 
specific angular momentum of the stars. 

We calculated Ar in our sample both radially (see Fig. as 
well as within one re (Arc, given in Tab. |^. We normalised the 
radial profiles to the bar radius to test the infiuence of the bar on the 
shape of the profile. Considering the small sample size and large 
variety of Hubble types and bar types, it is not surprising that we 
recover a variety of profiles. Yet more than 70% show a dip in Ar 
at around 0.2zb0.1 Rbar- The only galaxy which shows a clear offset 
of this dip is NGC 4262, where the stellar and gas velocity fields are 
clearly misaligned. 

This feature appears to be related to the double-hump in the 
velocity profile ( [Bureau & Athanassoula|[2005] ), in combination 
with the rise in a after the a—drop. Three galaxies (NGC 1015, 
NGC 2893, and NGC 5350) do not exhibit a double-hump nor 
a strong a—drop. In other galaxies (NGC 2712, NGC 4267, 
NGC 4262, NGC 5375), the hump feature in the stellar velocity 
alone seems to be strong enough to produce the drop in the Ar 
profile. In other cases, in particular NGC 4245, NGC 3485 and 
NGC 5701, the peak of the hump in the velocity profile coincides 
with the peak of the velocity dispersion profile (after the central 
drop). Thus, despite the distinct morphologies and inclinations, we 
observe a common behaviour and influence on the stellar kinemat¬ 
ics. This is likely produced by the bar, since these features seem to 
accumulate around a similar radius related to the bar length. This 
feature could probably be associated to the inner Lindblad reso¬ 
nance (ILR) (e.g. |Elmegreen|1994||Pfenniger & Norman|1990| l. In 
fact, in several of our galaxies, nuclear rings have been detected 
at those locations: NGC 2859 ( [Erwin & Sparke|2002j ), NGC 3504 
([Buta & Crocker[[1993| [Elmegreen et al.[[1997[l, NGC 4262 and 
NGC 4245 ([Comeron et aL|20Tolf 

We also tested the location of the dip as a function of the ef¬ 
fective radius of the bulge and the disc scale-length (derived from 
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Figure 10. Radial profiles of Ar for all galaxies, normalised to the radius 
of the bar. A rise and consecutive dip is observed at a similar position in the 
profiles of the majority of galaxies, around 0 . 2 ± 0.1 Rbar- 


Figure 11. Kinematic gas torque versus kinematic stellar torque. The dotted 
line indicates a 2.5-correlation. Grey galaxies show obvious offsets between 
gas and stellar kinematics. 


Table 3. Summary of the values obtained for A^e and the bar strength 
measurements: (1) NGC number, (2) (3) photometric torque, (4) error 

in the photometric torque, (5) kinematic torque, ( 6 ) error in the kinematic 
torque, (7) kinematic torque of the gas component, ( 8 ) error in the kine¬ 
matic torque of the gas component. As the error of the measure for A^e is 
negligible, we do not list it here. 


Galaxy 

(1) 

Arg 

(2) 

Qb 

(3) 

AQb 

(4) 

Qkin 

(5) 

AQkin 

(6) 

Qkin ,g 
(7) 

AQkin 

(8) 

NGC 1015 

0.25 

0.26 

0.074 

- 

- 

- 

- 

NGC 2543 

0.62 

0.36 

0.070 

0.39 

0.069 

1.1 

0.23 

NGC 2712 

0.65 

0.28 

0.044 

0.37 

0.044 

0.46 

0.07 

NGC 2859 

0.37 

0.17 

0.025 

0.22 

0.044 

0.94 

0.40 

NGC 2893 

0.06 

0.16 

0.020 

- 

- 

- 

- 

NGC 2962 

0.44 

0.14 

0.024 

0.080 

0.040 

0.62 

0.44 

NGC 3485 

0.52 

0.38 

0.064 

0.33 

0.084 

0.84 

0.43 

NGC 3504 

0.29 

0.26 

0.044 

0.24 

0.082 

1.5 

0.99 

NGC 4245 

0.33 

0.18 

0.020 

0.10 

0.033 

0.20 

0.19 

NGC 4262 

0.33 

0.07 

0.012 

0.14 

0.048 

1.8 

0.25 

NGC 4267 

0.24 

0.04 

0.013 

- 

- 

- 

- 

NGC 4394 

0.46 

0.23 

0.036 

0.23 

0.12 

0.52 

0.12 

NGC 4643 

0.28 

0.28 

0.069 

- 

- 

- 

- 

NGC 5350 

0.62 

0.44 

0.076 

- 

- 

- 

- 

NGC 5375 

0.47 

0.23 

0.044 

- 

- 

- 

- 

NGC 5701 

0.20 

0.18 

0.022 

0.08 

0.076 

0.42 

0.66 


the S^G) but did not find any correlation. This supports our sus¬ 
picion that this feature is related to the bar. Similar studies have 
already related outer ring radii to the bar sizes (e.g., [Perez, AguerrT| 
|& Mendez-Abreu||2012| | based on earlier studies and simulations 
(e.g. [Byrd et al.|199l]|Buta et al.|1995j l. |Comer6n et al.| ( [20T0] ) es¬ 
timate that the maximal possible extension of a nuclear ring should 
be located at 0.25 bar lengths. 


5 CONTRASTING BAR STRENGTH MEASUREMENTS 

In this section we compare bar strength measurements determined 
according to the descriptions given in Section [33] For the case of 
Qkin we additionally measure it for both the stellar and ionised-gas 
components. This comparison will allow us to establish the ruler 
that will be used in Sect.j^to evaluate the impact of bars of different 
strengths on different kinematic properties of our galaxies. 


5.1 Qkin from stars and ionised gas 

We calculated the kinematic torque (Qkin) from both our stellar and 
ionised-gas velocity maps (values given in Tab.[^. The comparison 
is shown in Fig. m The gas is expected to respond more strongly 
to the bar than the stars. Overall the torque derived from the gas 
velocity fields is significantly higher than the stellar one. For al¬ 
most all galaxies it is about 2.5 times greater than the value derived 
from the stellar velocity maps (dotted line). The grey points mark 
significantly higher values. Inspecting their gas velocity maps, we 
detect clear differences from the stellar velocity fields. In partic¬ 
ular NGC 4262 shows the highest value of Qkin,gas- The gas ve¬ 
locity field is counter-rotating with respect to the stellar velocity 
field. Therefore, a significant impact due to another process (e.g. 
galaxy interaction or close encounter) might be at work in this 
galaxy (see jVollmer, Huchtmeier & van Drielj2005| l. In NGC 3504, 
the gas velocity field shows extreme enhancements in the central 
regions which are not present in the stellar velocity field. It is not 
clear at this stage what is causing this difference. Overall, we find 
that the stars seem to be more stable and therefore the stellar kine¬ 
matic torque agrees better with the photometric torque (see § |5.2j ), 
whereas the gas is more susceptible to other processes, leading to a 
larger number of outliers. 
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5.2 Kinematic vs photometric bar strengths 


Figure[^(top panel) compares the kinematic (Qkin) versus the pho¬ 
tometric (Q^) torque measurements for the subsample of 10 galax¬ 
ies, where the kinematic method was possible. Despite large uncer¬ 
tainties, the correlation between the two parameters is obvious. This 
is confirmed by a measured linear Pearson correlation coefficient of 
0.83. The distribution of existing bar strengths within our limited 
sample is representative of larger samples of nearby galaxies ( |Lau-| 
[rikainen et al.|2004] ). Overall, early-type galaxies have lower values 
than the later types, confirming earlier results by e.g. |Buta et aT] 
( |2005| ); [Laurikainen et al.| < [2007] ). This could be a result of different 
factors: i) the infiuence of stronger spiral arms that still alter the mo¬ 
tions within the bar region (although we tried to avoid them in our 
analysis), ii) the presence of more gas in later types which is more 
responsive to the bar could also infiuence the stellar motions, iii) 
discs in earlier types are simply hotter leading to more random ver¬ 
sus ordered rotational motion. The dotted line in the figure indicates 
a one-to-one correlation. 


To investigate further the relation between stellar Qkin and 
, we have produced an extensive set of numerical simulations of 
barred galaxies following those in[Martinez-Valpuesta, Shlosman &| 
|Heller| { [2()06| ) and [Martinez- Valpuesta & Gerhard ( 2011| l. Here we 
use four simulations series, h, I 2 , h and I a, each one with a differ¬ 
ent disc-to-total ratio: 0.92, 0.62, 0.43, 0.29 respectively. This setup 
allows us to explore the effect of distinct dark matter haloes on the 
torque parameters. We analysed 1800 snapshots taken at different 
points in time of the bar evolution. In addition, we also varied the 
inclination and the position angle of the bar relative to the galaxy’s 
position angle to have different viewing angles and thus assess the 
infiuence of these parameters (see Appendix [B] for more details). 
The bar strength measurements of the simulations are presented in 
Fig. [T^ (bottom panel). 


We analysed the simulations in the same way as the obser¬ 
vations. We calculated from their simulated intensity distribu¬ 
tion and Qkin from their associated stellar velocity maps (see Ap¬ 
pendix]^ for details). Due to the lack of rK 2 o (k-band photomet¬ 
ric parameter) to infer the scale height hz (for the calculation of 
Q^) (jSpeltincx, Laurikainen & Salo|2()08] l, we applied the jde Grij^ 
( j 199 ^ relation for intermediate type galaxies which links the scale- 
height to the scale-length, assuming an exponential disc without 
truncations. The overall trend found is consistent with what we find 
with the observations. Nonetheless, distinct simulation series be¬ 
have systematically differently; the figure shows that higher disc 
fractions consistently lead to lower bar strengths, both in Q^^ and 
Qkin . Each simulation series exhibits low bar strengths, which cor¬ 
respond to snapshots in very early times in the bar formation. While 
1 1 soon seems to saturate and cannot grow stronger bars, the others 
do and saturate at later stages such that the strongest bars are found 
in the simulation series I a, the one with the highest dark matter 
fraction. 
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Figure 12. {Toppanel) Kinematic versus photometric torque (Qkin vs Q^^). 
The dotted line indicates a one-to-one correlation. Symbols are colour- 
coded according to Hubble type. The relation has a linear Pearson correla¬ 
tion coefficient of R=0.83. (Bottom panel) Comparison of Qkin and Q^ for 
the observations (black) and the four sets of numerical simulations {Ii, I 2 , 
1 3 and 74 ). We only plot stages of the simulations that exclude the buckling 
phases. 


6 THE EFFECT OF BAR STRENGTH ON GALAXY 
PROPERTIES 

In this section we try to understand if stronger bars affect the prop¬ 
erties of the host galaxy in a systematic way, focusing on whether it 
leads to stronger or weaker kinematic features. 


Given the good agreement between Q^^ and Qkin for our subset 
of galaxies along with the large number of simulations, we will use 
the photometric values determined from the S^G images for the bar 
strength values, because these are available for our entire sample. 
Our study also serves for verifying the technique and results of Q^. 


6.1 Relation with Hubble type 

Figure [T^ illustrates the already observed trend of Q^ with Hub¬ 
ble type (e.g., [Laurikainen et al.|2007| t, resulting in a linear Pearson 
correlation coefficient of 7^ = 0.96 for our sample (averaged val- 
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Figure 13. Bar strength as a function Hubble type. We represent our sam¬ 
ple averaged per Hubble type in comparison with the sample analysed in 
[Laurikainen et al.||2007) . 
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ues per Hubble type bin). In comparison with the observed trend 
found in [Laurikainen et aL] ( |2007| ), we cover slightly stronger bars 
throughout but conserving the trend, making our sample represen¬ 
tative of rather strongly barred galaxies. For a wider study of as 
a function of Hubble type based on the S^G sample, please refer to 
Dfaz-Garcfa (2015, submitted). 

In Fig. [M we depict the position of the dip in the Ar profile, 
depicted in Fig as a function of bar strength. The dip is not ob¬ 
served in all galaxies, therefore we only show those, which exhibit 
this feature. Is this dip feature related to inner structures such as 
nuclear rings? We plot in gray positions of rings from the AINUR 
sample ( [Comeron et al.|201^ , including three galaxies (NGC 2859, 
NGC 3504 and NGC 4245) that we share. Evidently, ring and dip 
positions do not correlate and we do not find any mathematically 
significant trend for either sample. 

Fig. shows again the position of the dip in the Ar pro¬ 
file, now as a function of light concentration R90/R50. Apart from 
NGC 4262 - already found not to follow other observed trends, 
probably due to a recent interaction - the galaxies seem to follow 
a downward trend: the more concentrated the bulge, the closer is 
the dip feature towards the center. This could be directly related to 
the bulge: in our simulations with more concentrated bulges, we 
also find that the ILRs are located closer to the center. It could also 
mean that these features are more evolved in time, supporting the 
prediction of the migration of nuclear rings towards the center (e.g.. 


Knapen et al.|1995[|Fukuda, Habe & Wada|2000[ [Regan & Teuben 


2003 [van de Ven & Chang|2009| >, also recently observed by Pinol 
Ferrer et al.[ ( [2014[ l. 

This trend is only mildly observed for the values of the AINUR 
sample, taking their measured ring radii as a comparison, because 
no Ar profiles are available for that data. We further determined 
the position of ilLR and oILR by a simple linear approximation 
analysing U-curves and estimates for the bar pattern speed. Neither 
the position of the dip in the Ar profile, nor the ring radius are found 
at the exact same position as these resonances (in a forthcoming 
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Figure 14. Position of the dip in the Ar profile as a function of bar strength 
based on the profiles obtained in Fig[T^ Grey points indicate ring positions 
(as bar length fractions) also as functions of the bar strength, both as mea¬ 
sured in the AINUR sample (Comeron et al.[2010) . 
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Figure 15. Position of the dip in the Ar profile as a function of light con¬ 
centration R 90 /R 50 based on the profiles obtained in Fig[^ Grey points 
indicate ring positions (as bar length fractions) also as functions of the con¬ 
centration, derived in the same way as for our sample. 
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paper, we will verify this by a more robust calculation of bar pattern 
speeds). 


6.2 Influence of the bar strength on the global position angle 

To quantify the influence of the bar on the global velocity map, 
we analysed the difference between the photometric and kinematic 
position angles as a function of the bar strength. We compared the 
photometric PA with the stellar and gas kinematic PA, as well as 
the difference between the stellar and gaseous components. This is 
shown in Fig.[^ We measured the kinematic PAs directly from the 
velocity maps following [Barrera-Ballesteros et al.[ ( [20T4] ). 
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Figure 16. Photometric and kinematic position angle differences between 
stars and ionised gas as a function of bar strength. Top: difference between 
global photometric and stellar kinematic PA. Middle: difference of photo¬ 
metric to kinematic gas PA. Bottom: difference between the stellar kinematic 
and gas kinematic PA. Representative error bars are indicated in the top right 
corners. 


No trends are observed with bar strength in any of the three 
cases. In previous studies including barred galaxies (e.g. , [Falcon 


Barroso et al.|2006[|Fathi et al.|2009[[&ajnovic et al.|201 l[|Barrera- 


Ballesteros et al.|2014 l, the observed misalignments were neither 

strong nor systematic. Misalignments were found mainly in systems 
with complex kinematics (non-regular rotators), systems in interac¬ 
tion and only in some barred galaxies, but their amplitudes largely 
depended on the FoV. The detailed study by |B arrera-B allestero^ 
jet al.j ( |2014j ) further concludes that morphological substructures 
only influence the redistribution of angular momentum, but the 
global kinematics such as rotation are driven by the overall disc 
mass. Here, we also only detected one galaxy, NGC 3504, with a 
larger difference between the photometric and kinematic PAs. The 
ionised gas is known to react strongly to the bar producing a twist 
in the zero-velocity curve (e.g., [Peterson & Huntl^|1980| [Em-[ 
[sellem et al.[2006[ l. However, for the gas too, we only detected large 
misalignments in NGC 4262. The overall misalignment is slightly 
larger than that for the stellar kinematic/photometric PA difference, 
but not significantly. Finally, the difference between the two kine¬ 
matic PAs (stellar and ionised gas PA) is found to be equally small 
and not correlated with the bar strength, confirming the results from 
earlier studies (as described above). 


6.3 Stellar kinematic features related to bar strengths 

Bars can be depicted as engines that, on the one hand, drive gas 
towards the central regions and consequently nourish star forma- 
tion (e.g.[Elmegreen & Elmegreen[1985][1989[[Erwin[2005][Ellison[ 
[et al.|201l|l and, on the oth er hand, support radial motions of stars 

(e.g., [Minchev & Eamaey[[2010[ [Brunetti, Chiappini & Pfenniger| 

[201 1[ |. Due to these factors, they are natural triggers of changes in 
the centre of galaxies, and in the stellar velocity dispersion in partic¬ 
ular. We do not find, however, any trend between the central velocity 
dispersion and the bar strength. This reinforces the picture that the 
central stellar velocity dispersion is determined by global galaxy 
properties. At least, it does not vary significantly due to morpho¬ 
logical substructures in a systematic way, except for the occasional 
central cr—drop. 

While simulations predict and find a significant influence of 
the bar on the host galaxy in various ways (e.g., [Martinez- Valpuesta, 
Shlosman & Heller[[2006[ [Athanassoula, Machado & Rodionov 
201 3[ [Sellwood[[2014[ l, we only find mild signatures on the kine¬ 
matic maps in our sample, such as the proposed double-hump ro¬ 
tation curve and occasional a—drops. Despite this lack of major, 
bar-induced alterations in global galaxy kinematic parameters, we 
detect some relation between those subtle kinematic features and 
the bar strength. It is thus logical to assume a connection between 
those features and the bars. 

Since double-humps and cr—drops exist commonly among 
barred galaxies, we tested their amplitude in relation to the bar 
strength. In other words, would stronger bars produce stronger 
humps or deeper drops? We quantified the strength of the hump 
by the difference of its inner peak and consecutive drop, calling this 
parameter AV. We further normalised this value by the maximum 
rotation - corrected for inclination - that we could detect for each 
galaxy. We are aware that asymmetric drift could attenuate this sig¬ 
nal slightly, but do not expect a major change for the trend observed. 

For the velocity dispersion we determined the amplitude of 
the central cr—drop if present, again normalised by the maximum 
velocity dispersion (following [Peletier et al.[[2012] l. We chose to 
compute these quantities at the position angle where the signal was 
stronger. Since the hump in the velocity is seen strongest along the 
major axis, we took the profile along that axis. We chose to take the 
velocity dispersion profiles along the bar major axis, because the 
drop is most pronounced along that direction. 

Figure [T^ shows the results. We identified galaxies with very 
low inclinations in light grey, galaxies with intermediate but still 
low inclinations or larger uncertainties in their velocity fields in 
grey, and reliable points in black. As our sample is very small, we 
did not discard any points, but indicate that we are conscious about 
the bias introduced by measuring at different inclinations. As incli¬ 
nation effects in the velocity dispersion are very difficult to char¬ 
acterise (i.e., it depends on the projection of the velocity ellipsoid 
being probed and anisotropy), we did not attempt any correction. 

We find tentative evidence that stronger bars produce stronger 
humps in the velocity profile. After the inclination correction, the 
low-inclination galaxies also follow this trend and we obtain a linear 
Pearson correlation coefficient of R = 0.57. Discarding unreliable 
galaxies (light gray points), the correlation coefficient increases to 
R = 0.76. As the hump could sometimes be distinguished better in 
the ionised gas, we also determined these parameter in the gas ve¬ 
locity profile (not shown here). The results follow the same trend. In 
the bottom row of Fig. EZl we show the measurement of the mag¬ 
nitude of the cr—drop. Stronger bars produce a stronger cr—drop 
features. We obtained an overall linear Pearson correlation coeffi- 
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Figure 17. Different parameters detected in the kinematic maps of barred 
galaxies as a function of bar strength. Top: Magnitude of the velocity hump 
along the major axis rotation curve, normalised by the maximum rotation 
velocity and corrected for inclination (see text for details). Bottom: Magni¬ 
tude of the cr—drop, normalised by the maximum dispersion (i.e. difference 
between central drop and highest surrounding velocity dispersion). In both 
panels, the dashed line indicates a linear fit to the black and dark gray points. 



cient of R = 0.74. All except one of the galaxies without a central 
dispersion drop have values below 0.15. The lack of this drop 
feature seems to be most evident in galaxies with weaker bars. 


6.4 Stellar angular momentum as a function of bar strength 

We now inspect the influence of the bar on the integrated angu¬ 
lar momentum within one effective radii (Arg). This is shown in 
Fig. The values obtained for Arg are consistent with the values 
found in the literature for barred galaxies (e.g. [Bender, Saglia &| 
|Gerhard|1994[|Krajnovic et al.|2008||201 1 j i. These studies, however, 
do not include bar strength measurements. We observe an increas¬ 
ing value of Arg with bar strength. This is somewhat connected 
to Hubble type, because the later-types in our sample display the 



Qb 


Figure 18. Apparent stellar angular momentum within one effective radius 
i jEmsellem et al.|2()07) as a function of bar strength. The different Hubble 
types are shown in different colours. 


largest Arg values. The high Arg values observed in the later-type 
galaxies are likely due to the higher fraction of the disc, and thus 
high rotation, included within the one effective radius aperture. 

As there is angular momentum transfer between the bar, 
disc and outer halo (e.g., [Combes & Elmegre^|1993[ [Martinez^ 
|Valpuesta, Shlosman & Heller [[2006[ l, the value of Arg should be 

higher in barred galaxies compared to their non-barred counter¬ 
parts. [Athanassoula[ ( [2003[ ) showed, however, that while angular 
momentum is transferred to the disc, the bar also slows down and 
therefore contributes to a decrease in Arg. The current available 
Arg values in the ATLAS3D ( [Emsellem et al.[[2011[ l or CALIFA 
samples ( [Falcon-Barroso et al.[2014[ ) do not indicate distinct values 
for barred and non-barred galaxies. 


6.5 Bars as drivers of radial motions 

Bars have been studied as a major driver of radial mixing for a long 
time (e.g.[Friedli & Benz|199^ , but spiral arms (e.g., [Sellwood &[ 
Binney|2002[l or the comb in^on of their resonance overlap (e.g., 
Minchev & Famaey|2010[[Shevchenko|2011[ l are also held respon¬ 
sible for an increase of this. Investigating the latter, [Brunetti, Chia^ 
[pini & Pfenniger[ ( f2011[ ) found that kinematically hot discs are not 
as efficient environments and exhibit less radial mixing than kine¬ 
matically colder ones. We investigate the level of radial motions in¬ 
duced by bars in our sample by studying the radial gradients of the 
stellar velocity dispersion. The expectation is that barred galaxies 
would display flatter gradients than those measured in non-barred 
systems. 

We start by comparing velocity dispersion profiles along the 
major and minor axes of the bar. In earlier literature studies, ve¬ 
locity dispersion profiles were investigated typically only along the 
major axis of the bar (e.g. [Perez, Sanchez-Blazquez & Zurita|2009[ l. 
For our sample, this is shown in Fig.[19|for one galaxy as an exam- 
ple (similar plots for other galaxies are presented in Appendix [C]). 
We show the overall stellar velocity dispersion map for reference 
and, next to it, the radial profiles along the major and minor axes 
(minor axis radii are corrected for inclination) of the bar extracted 
from this map. The profiles along the axes are generally overlap¬ 
ping, hence we do not observe any increase along the bar major 
axis. In one third of our sample, however, we observe a mild dif- 
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Figure 19. Stellar velocity dispersion map of NGC 4643 and associated profiles. {Left) Stellar velocity dispersion map with bar major and minor axes indicated 
with dashed lines. {Middle) Averaged profiles, corrected for inclination. Dashed line indicates the bar length. {Right) Averaged profiles along isophotes (first 
point corresponds to first isophote seen on the top left image). 


ference around the central parts (also seen in the example), with 
the major axis showing a higher dispersion. This is probably linked 
to the aforementioned kinematic substructures such as inner discs 
or rings, possibly a result of barred secular evolution. The major 
axis profiles observed in |Pdrez, Sanchez-Blazquez & Zurita| ( |2009| ) 
show a similar behavior to ours, but no minor axis measurements 
have been performed in that work. 

As bars are structures seen in the photometry, we decided 
to also trace the profile comparing their points along the same 
isophotes. The isophotal profiles reveal a larger overall a along of 
the major axis than the minor axis, and not just in the central parts. 
The fact that the velocity dispersion further out is higher along the 
minor axis, compared with the same isophote on the major axis, 
shows that the dispersion of the bulge - traced by the minor axis 
- dominates strongly, regardless of the prominent bar seen in the 
photometry. It shows nonetheless that the kinematics of the bar is 
significantly different than the bulge and it is more similar to the 
disc, because at the outer end of the profiles, reaching the disc, val¬ 
ues along the major and minor axes start coinciding again. 

Figure shows a relation between the outer gradient of the 
stellar velocity dispersion and the bar strength. The gradient is shal¬ 
lower for stronger bars. Nevertheless, the trend is based on only 
very few points, in particular in the low bar strength regime. The 
presence of a bar can cause enhanced radial motions which perturb 
the system. Thus the orbital mixing increases which in turn can lead 
to higher dispersion and shallower gradients. This could be the rea¬ 
son for the observed flattening of the gradients with higher . Ad¬ 
ditionally, we measured the gradient not only along the major axis 
of the galaxy (black points), but also along the minor axis (blue 
points). If the bar would significantly flatten the gradient along the 
major axis, the minor axis values would be expected to show steeper 
gradients. The results, however, show a scatter of shallower as well 
as steeper gradients along the minor axes compared to the points 
measured along the major axes and we cannot identify a systematic 
behavior. 

The measurement of the maximal radial velocities that we use 
to calculate the kinematic torques also indicate the average radial 
displacement. As mentioned in previous sections, the average radial 
motion is much higher in the gas than for the stars. Stars move at 
velocities between 10 and 60 km/s radially, corresponding to about 
10-60 kpc/Gyr or 0.09 - 0.88 when normalised with the rotation at 
~ re , a value which is similar in magnitudes to what we And with 
our own simulations (e.g. in Ia it is around 0.2) whereas the gas 
move at 40 - 100 km/s and in extreme cases such as NGC 4262 at 
more than 300 km/s. The latter is most probably due to an outer 



Figure 20. Gradient of the outer velocity dispersion, ignoring the central 
regions, as a function of bar strength. Blue points indicate profiles taken 
along the bar minor axis, while black points show the ones along the bar 
major axis. 

influence such as an interaction (Vollmer, Huchtmeier & van Driel| 
|2005| ). Nevertheless, the stronger effect on the gas than on the stars 
has been seen in numerous simulations ( [Athanassoula, Machado] 
& Rodionov|2013[| Kubryk, Prantzos & Athanassoula|2013| e.g.). 

Kubryk, Prantzos & Athanassoula| ( |2014] | And a particular influ¬ 
ence of the bar-induced radial inflow on the gaseous profile. Fur¬ 
thermore, |Maciejewski, Emsellem & Krajupyid] ^20l2} obtain val¬ 
ues which are in the range of the ones we recover using the same 
method. The recent work of |Goz et ^ ( |2015| ) analyzing two simula¬ 
tions of barred galaxies resulting from N-body-fSPH cosmological 
simulations, shows one case with significantly higher radial mo¬ 
tions (around 150 km/s) whereas in the other case the magnitude is 
comparable to what we observe (around 30 km/s). 

In a large number of simulations (e.g., |Minchev et al.||2012| | 
bars are found to be the most efficient driver of radial migration, in 
particular through their corotation resonance. We have not yet de¬ 
termined the radius of corotation for our sample, but plan to com- 
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pute this in a forthcoming paper. This parameter, together with the 
stellar populations will further complete the picture of bar-induced 
mixing. In particular in Paper II of this series, we will assess the 
impact of the radial motions determined here, on the stellar popula¬ 
tion properties, which will allow us to shed light onto radial migra- 
tio n effects (e.g.,|Friedli, Benz & Kennicutt|1994l Haywood|2008[ 


[Roediger et al.|2012[|Kubryk, Prantzos & Athanassoula|2013| l. 


7 SUMMARY AND CONCLUSIONS 

We present the BaLROG sample of 16 barred galaxies of different 
Hubble types, spanning the typical bar strengths found in the local 
Universe. Our large mosaics with the integral field unit SAURON 
cover the bars out to the radius where the disc begins to dominate, 
at a spatial resolution of typically 100 pc. For every galaxy we also 
use Spitzer observations from the S^G survey of nearby galaxies 
( |Sheth et al.|2010| ) to determine several photometric parameters, as 
well as to derive the bar strength . From the velocity maps, we 
calculate radial and tangential velocities to compute the bar strength 
based on the kinematics, Qkin- Our aim is to establish a reliable 
yardstick, namely bar strength, to probe the influence of the bars on 
different parameters of the host galaxies. 

In this paper we focus on the kinematics of the galaxies, de¬ 
riving stellar and gas velocities and velocity dispersions, ha and 
h4 Gauss-Hermite moments and the stellar angular momentum Ar 
and carefully comparing to a large set of N-body simulations. The 
analysis of our observations leads to the following results and con¬ 
clusions: 


• Bars do not strongly infiuence the global kinematics of their 
host galaxies, regardless of their strength. Our work confirms pre- 


vious studies (e.g. [Falcon-Barroso et al.||2006 

[Fathi et al.||2009| 

[Krajnovic et al.||2011| Barrera-Ballesteros et a 

l.||2014|l and shows 


the lack of strong kinematic misalignments between the galaxies’ 
photometric and kinematic axes. 


• Bars do have an infiuence on more subtle kinematic features, 
especially in the inner regions of galaxies. We detect double-hump 
velocity profiles and velocity dispersion drops (e.g., [Bureau &| 
|Athanassoula|2005] l, which increase in intensity with increasing bar 
strength. 

• We find evidence for the presence of inner structures such as 
inner rings or discs in about 50% of our sample. These features are 
detected from the anti-correlation between hs and and V/cr within 
the effective radius of the galaxies (?^0.1 bar lengths). 

• The derived Ar profiles show a dip at 0.2±0.1 Rbar, which we 
suggest is connected to the presence of inner substructures. 

• We also derived the integrated angular momentum within one 
effective radius (Arg) and find that galaxies with stronger bars ex¬ 
hibit a higher Arg value. This may be a secondary effect of late-type 
galaxies, because they are more rotationally supported and thus also 
host stronger bars. 

• We developed a new method to determine the bar strength 
from stellar or ionised gas velocity maps (Qkin). This method re¬ 
lies on the extraction of the ratio of radial and tangential velocities 
using the technique developed by [Maciejewski, Emsellem & Kra-| 
[jnovic| ( [?012j ). Values of this parameter agree well with independent 
measurements obtained from imaging, , e.g.,|Laurikainen & Saloj 
( |2002| ), and predictions from numerical N-body simulations. 

• Bar strength values measured from ionised-gas kinematic 
maps are a factor ~2.5 larger than those determined from the stellar 
kinematic maps. 


• We observe a flattening of the outer stellar velocity dispersion 
profiles with increasing bar strength. 

These results suggest a complex influence of bars in nearby 
galaxies, especially affecting central regions. We do not observe a 
significant influence on global properties, but bars seem to affect 
only on small scales. The gas is clearly more strongly affected, re¬ 
flected in higher gaseous than stellar torques. In our sample we de¬ 
tect a difference between bars in early and late-type galaxies hinting 
towards a different mechanism, maybe due to the presence of higher 
and lower gas fractions. To better answer these questions and deter¬ 
mine time scales, we will investigate the stellar populations of these 
galaxies in detail in BaLROG II (Seidel et al., in prep) and also de¬ 
termine their corotation radii and pattern speeds in relation with 
their dark matter fractions (BaLROG III). 
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APPENDIX A: DESCRIPTION OF THE SIMULATIONS 

In the following we illustrate the give a more detailed description 
of the simulations and how we used to obtain the torque measures. 
Eigure [AT] shows as an example the intensity and velocity maps of 
one of the snapshots of the /s simulation at an intermediate time 
step, at an inclination of 30 degree and a bar-to-line-of-nodes-PA of 
55 degree. EigurejA^ shows the corresponding obtained radial and 
tangential velocities. In Eig. |A3| the full time series (all 1800 snap¬ 
shots) of this particular simulation is shown. The early very strong 
peak is associated with the first buckling of this bar (not stable in 
time). The shape of the obtained curve indicating the kinematical 
torque values resembles the curve of the A2 value measured dur¬ 
ing the simulation and also gives an indication on the bar strength. 
Eor the other simulations series, these measurements along with the 
shape of the final curve differs significantly. 

The variations of Qkin with inclination and PA (shown with 
different colors and symbols respectively in Eig. |A3| are due to the 
fact that the assumptions of a thin disc and stable bar are not 100% 
correct. This is the case especially during the buckling event of the 
bar and during the later stages of the evolution, which is to be ex¬ 
pected due to the thickening of the bar. If it was a perfect measure¬ 
ment, all points should overlap vertically, since we simply rotate 
the simulated galaxy slightly in order to achieve its different orien¬ 
tation. Since the buckling event is a short moment during bar evolu¬ 
tion, we can safely exclude these points from our comparison as the 
likelihood to find a bar in the buckling phase is rather low. Without 
these points and below t=200, the measurements agree rather well. 
A more extensive test can also be found in the Appendix [B] 

To perform the analysis for the comparison of the two torque 
measurements, we calculated the values of Qb (in addition to those 
of Qkin) for all sets of simulations in the same way as done for 
the observations. In Pig. we presented the final results showing 
the observational measurements overlaid on the four simulations, 
A, h, h and h (level of disc-to-total ratios: 0.92, 0.62, 0.43, 0.29 
respectively). 

As seen in Pig.[^ none of the simulation series coincides per¬ 
fectly with all the data points, but instead form a continuum. It is 
reassuring however that the simulation series fall onto the same re¬ 
lation found as for the observations. 

Simulation h (92% disk, shown by the green points) best 
represents the lower bar strengths. With decreasing disc content, 
the both the photometric and kinematic bar strengths increase. 
We also distinguish a particular behavior for I 4 : compared to the 
other sets, this series seems to adapt to all ranges starting from 
very weak to very strong bars. This is the simulation with the 
lowest disc percentage within the bar region initially and also after 
halo relaxation and is at every point in time submaximal. Overall, 
we find a trend of stronger bar development with increasing dark 
matter halo fraction in our simulations. 

Comparing the velocity maps of the three different simula¬ 
tions with the observations, we find good agreement, especially in 
the case of Ii . Simulations I3 and I4 in particular develop strong 
distortions in the stellar velocity fields which are not as pronounced 
in our observations. However, the bar strength of A never reaches 
a higher value than about 0.2-0.3, unlike the observations. To reach 
higher values, we needed to increase the halo fraction. 
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Figure A3. The kinematical torque evolution for one complete simulation, 
here /a as an example, indicating different inclinations and PAs. The dotted 
line indicates the point in time where the bar has reached its full strength 
and further evolution is not reliable, hence points to the right of it will not 
be considered in our analyses. 


Figure Al. Intensity (in logarithmic units) and velocity maps (translated in 
km/s) of one of the snapshots of the /a simulation at an intermediate point of 
the bar evolution t=150 at PA = 35 and inclination = 30. The bottom panels 
show the prepared maps which can be symmetrised by the code, with the bar 
in a horizontal position and adjusted field extensions. In those, we overlay 
the isophotes. 




Figure A2. The obtained radial and tangential velocities for the simulated 
galaxy above, but now symmetrised. 


higher values in all three parameters, whereas for higher PAs, high 
inclinations result in higher values overall. The spread here is less. 

For the other two simulations with higher DM content (/s and 
I a), the effect of inclination is always the same despite the distinct 
PAs: a lower inclination results in higher values in the three mea¬ 
sured parameters. Again, the spread is less at higher PAs (except for 

Qkin in F4). 

The fact that the influences are similar in spite of the different 
measurement methods probably helps to produce the observed re¬ 
lation between them. It is important to bear the influence of these 
parameters in mind when checking the observations: low inclina¬ 
tions might lead to higher values, and, depending on the PA and 
DM fraction, high inclinations can also lead to lower values. 

We also compared Qb with the A2 values directly from the 
simulations and And that Qb resembles A2 very well. 

Overall, values of Qb are expected to be higher at lower in¬ 
clinations whereas values of Qkin should be higher at higher incli¬ 
nations as motions can be better measured with increased inclina¬ 
tion. Our tests, however, suggest that in almost all cases, the chosen 
methods indicate a lower limit for the bar strength, in particular in 
the case of Qb. 


APPENDIX B: INFLUENCE OF THE INCLINATION AND 
PA ON THE BAR STRENGTH MEASUREMENTS 

We used four simulation sets (/i, I 2 , h and I a) to test the influence 
of PA and inclination on the bar strength measurements, namely 
Qb, Qkin and A2, on a large enough sample. Figure [ bT] illustrates 
those tests with the example of the I 2 simulation. 

In the majority of cases the influence of these two parameters 
causes a consistent change in any of the different bar strength mea¬ 
surement methods. The PA influence always shows a clear trend: 
larger PAs result in higher strength values overall. Furthermore it 
causes less spread within the distinct inclinations, especially for Qb 
and A2. 

The effect of inclination is two-fold: in the simulations with 
higher disc-percentages (/i and I 2 ), its effect is reversed for low and 
high PAs. For low PAs, we detect that a lower inclination results in 


APPENDIX C: COMPLETE KINEMATIC MAPS FOR 
STARS AND IONISED GAS 

We show maps of the stellar and ionised-gas kinematics for the en¬ 
tire BaLROG sample of galaxies in flgures [CT] to |C16| In each flg- 
ure we show different maps of each galaxy, top to bottom and left to 
right: first row: (i) S^G image of the galaxy with an estimate of the 
final SAURON mosaic and the number of pointings indicated in the 
left lower corner, (ii) fundamental parameters of the galaxy along 
with the systemic velocity, inclination, stellar angular momentum 
within one effective radius and the bar strengths measured; second 
row: (i) surface brightness derived from the SAURON cube (col¬ 
lapsed in wavelength, shown in logarithmic scale), (ii) stellar mean 
velocity V (in km s“^), (iii) stellar velocity dispersion a (in km 
s“^); third row: (i) flux of the ionised gas, based on [Olll] (shown 
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Figure Bl. For a simulation with 60% disc we plot Qb and Qkin versus time for different inclinations and PAs: left: PA=25, right: PA=55. Top: Qb and 
bottom: Qkin - The elevation at early times is due to the buckling event in the bar evolution. 


as square-root-scaled), (ii) mean radial ionised gas velocity, (iii) 
ionised gas velocity dispersion (in km s“^); fourth row: (i) Gauss- 
Hermite moments h3 and (ii) h4; fifth row: (i) major and minor axis 
rotation curves of the stellar velocity, (ii) radial profile (inclination 
corrected) of the stellar velocity dispersion of the major and minor 
axis of the bar, (iii) isphotal profile of the stellar velocity dispersion 
along bar major and minor axis. The cut levels are indicated in a 
box on the right-hand side of each map. 
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Figure C2. Figure [CT] continued. 
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Figure C3. Figure[cT]continued. 
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Figure C4. Figure [CT] continued. 
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Figure C5. Figure [CT] continued. 
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Figure C6. Figure [CT] continued. 
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Figure Cl. Figure [CT] continued. 
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Figure C8. Figure[CT] continued. 
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Figure C9. Figure[cT]continued. 
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Figure CIO. Figure [cT] continued. 
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Figure Cll. Figure[cT]continued. 
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Figure C12. Figure [cT] continued. 
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Figure C13. Figure [cT] continued. 
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Figure C14. Figure [cT] continued. 
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Figure CIS. Figure[CT]continued. 
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Figure C16. Figure[cT]continued. 
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